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Upon increasing the electron density in a quantum wire, the one-dimensional electron system 
undergoes a transition to a (jruasi-one-dimensional state. In the absence of interactions between 
electrons, this corresponds to filling up the second subband of transverse quantization, and there 
are two gapless excitation modes above the transition. On the other hand, strongly interacting one- 
dimensional electrons form a Wigner crystal, and the transition corresponds to it splitting into two 
chains (zigzag crystal). The two chains are locked, so their relative motion is gapped, and only one 
gapless mode remains. We study the evolution of the system as the interaction strength changes, 
and show that only one gapless mode exists near the transition at any interaction strength. 

PACS numbers: 71.10.Pm 



Transport properties of quantum wires have attracted 
much attention over recent years [H, 0, 0, 0, 0, 0, 0] ■ Due 
to their quasi-one-dimensional structure, conductance is 
expected to be quantized in units of the conductance 
quantum Go = 2e 2 /h, where e is the elementary charge 
and the factor of 2 accounts for spin degeneracy. This 
property of non-interacting electrons is insensitive to the 
inclusion of interactions within the Luttinger liquid de- 
scription. However, a number of experiments show devia- 
tions from perfect conductance quantization, such as the 
so-called 0.7-structure observed below the first plateau 

0, 0, 0, 0] ■ These experiments have stimulated much 
theoretical interest in the physics of one-dimensional con- 
ductors not captured by the Luttinger-liquid theory, such 
as that of the spin-incoherent regime characterized by 
very weak coupling of the electron spins [1, 0, [l(| • Here 
we consider another important problem in this category: 
the transition of the one-dimensional electron system in 
a quantum wire into a quasi-one-dimensional state. 

Whether or not an electron system can be viewed as 
one-dimensional crucially depends on the strength of in- 
teraction 111 ]. In the absence of interactions, 
occupy subbands of transverse quantization 
system is one-dimensional until the chemical potential 
reaches the bottom of the second subband. On the other 
hand, at strong interactions, the electrons form a Wigner 
crystal, and the subband picture is no longer applicable. 
The system remains one-dimensional until the interaction 
energy overcomes the confining potential, and the crystal 
splits into two chains, forming a zigzag structure |l2j. 

There is an obvious difference in the behavior of the 
system in the vicinity of the transition between the limit- 
ing cases of the Wigner crystal at strong interactions on 
the one hand and non- interacting electrons on the other. 
In a zigzag Wigner crystal, the two chains are locked, 
and only one gapless mode (the plasmon) remains. By 
contrast, in the case of non-interacting electrons the two 



electrons 
and the 



subbands are independent and therefore represent two 
gapless modes. In this paper, we address the fate of the 
gapped mode in the vicinity of the transition as the in- 
teraction strength varies. In particular, we show that a 
gap exists at any interaction strength, Fig. [TJ However, 
the nature of the transition changes as the interaction 
strength is varied. 

For simplicity we consider spinless electrons and as- 
sume that they interact via long-range Coulomb repul- 
sion, 
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Here are the two-dimensional position vectors of the 
electrons, and e is the dielectric constant of the material. 
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FIG. 1: Phase diagram of the electronic system in a quantum 
wire. The interaction strength can be tuned by varying the 
confining potential. In the vicinity of the transition to a quasi- 
one-dimensional state (the lower line) the system supports 
only one gapless excitation mode at any interaction strength. 
The state with two gapless excitation modes is expected only 
at relatively weak interactions. 
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If the electrons in the wire are confined to one di- 
mension by a strong external potential V CO nt(yk), their 
physics is controlled by the one-dimensional electron den- 
sity n e . Since at n e — ► the kinetic energy of an elec- 
tron ~ (ti 2 /m)n 2 e scales to zero faster than the inter- 
action energy ~ (e 2 /e)n e , at low densities the Coulomb 
repulsion dominates. In this limit electrons behave classi- 
cally. In order to minimize their mutual interaction, they 
form a periodic one-dimensional structure — the so-called 
Wigner crystal. At small but finite density, quantum fluc- 
tuations smear the long-range order [l3j |. but the short- 
range order remains as long as the distance between elec- 
trons n~ x is greater than the Bohr radius as = eh 2 /me 2 . 

The above picture is valid if the width w of the wire 
is small, w <C as- In wider wires, or, equivalently, at 
stronger interactions, the opposite regime w 3> o,b can 
be achieved. In this case the electrons may form a two- 
dimensional structure while remaining essentially classi- 
cal. The structure of the Wigner crystal in this regime 
can be studied in detail (cf. Ref. Il2l ). if the confining 
potential is quadratic, 
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Here the frequency Q determines the width of the wire, 
w ~ \J h/mQ. The positions of all electrons are found 
by minimizing the energy Vint + V con { over keeping the 
one-dimensional electron density n e fixed. The geometry 
of the classical crystal is controlled by the dimensionless 
electron density v = n e ro, where ro = (2e 2 /emQ 2 ) 1 / 3 is 
the sole length scale of the problem ■ At v below the 
critical value v c = {/4/7C(3) « 0.780 the energy is mini- 
mized when electrons form a one-dimensional crystal, in 
which Xk = k/n e and yk — 0. At v > v c the crystal 
splits into two rows. The distance between rows van- 
ishes at the transition; just above the critical value, when 
5v = v - v c -C 1, it grows as c = ( v /24/93C(5)/V 2 )v / 5zJ 
(in units of ro). 

Let us consider the low-energy phonons in the zigzag 
Wigner crystal. Regardless of the density, the crystal has 
the usual plasmon excitation with acoustic spectrum. In 
the limit of zero wavevector, this excitation corresponds 
to translation of the crystal along the wire, Sxk = r\, 
5yk = for any k. In addition, at the zigzag transition 
point a transverse soft mode appears, for which Sxk = 
and Syk = (— l) fc </?. One can easily show that near the 
zigzag transition, when Sis ^ 1, the coupling of the two 
low-energy excitation modes is weak, and they can be 
treated separately. The action describing the soft trans- 
verse mode takes the form 

5- = Ah. / -^ f dr dx \(d T ip) 2 + (d x ipf -8vip 2 + ip 4 ] . 

V a B J ' 

(3) 

Here x, r, and the field ip have been rescaled so as 
to yield the simplest form of the action possible. This 



form, as well as our results from this point on, are 
not sensitive to the exact shape of the confining poten- 
tial. In case of the parabolic potential © the constant 
A = [7C(3)] 3 / 2 \/In2/31C(5). 

Above the classical transition point, i.e., at bv > 0, the 
transverse mode becomes unstable. This corresponds to 
the formation of the zigzag structure. To stabilize the 
system we keep the quartic term (p 4 . Quantum fluc- 
tuations affect both the position and the nature of the 
phase transition. To determine their effect, it is helpful 
to fermionize the problem. To this end we rediscretize the 
coordinate x and consider a set of particles moving in a 
double- well potential — Xip 2 + ip 4 with nearest neighbor 
interaction (ipj — ifj + i) 2 between them. At sufficiently 
large A each particle is almost completely localized in 
one of the minima, and its position can be described by 
a spin operator, ipj = ±a/A/2 = ^X/2a Z y In terms of 
these pseudo-spin variables the Hamiltonian contains two 
terms: —t ctJ describing tunneling between the two 
minima of the double-well potential and —vJ2j a j a j+i 
accounting for the nearest-neighbor interactions. Rotat- 
ing a x — * — <i z and applying the Jordan- Wigner transfor- 
mation, one obtains the Hamiltonian 



Ui = [ 2ta ho - v ( a ] - a j) ( a ]+i + 
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The above procedure has enabled us to convert the 
bosonic problem (J3J) to that of non-interacting spinlcss 
fcrmions Q. Since the number of fermions is not con- 
served, the Hamiltonian should be diagonalized by per- 
forming a Bogoliubov transformation. As a result one 
easily finds that the excitation spectrum of the Hamil- 
tonian (U|) has a gap A that vanishes when t = v. We 
identify this point with the phase transition from the 
one-dimensional state of the wire at t > v when all the 
fermionic states in the Hamiltonian Q are empty, to the 
quasi-one-dimensional state at v > t, in which fermionic 
states describing the transverse degrees of freedom in the 
wire are filled, but possess a spectral gap. Near the tran- 
sition the gap behaves linearly, A = 2\v — t\. 

In experiments with quantum wires, the transition 
from a one-dimensional to a quasi-one-dimensional state 
is observed when the chemical potential \i of electrons is 
tuned by applying a voltage to the gate controlling the 
electron density. The parameters t and v of the Hamilto- 
nian (HJ) are expected to be non-singular functions of /i. 
The transition occurs at the critical value fi c , defined as 
a solution of the equation t(/i) = v(fi). The gap in the 
excitation spectrum is then linear in the distance from 
the transition, 



A oc \fji - n c 
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To better understand the nature of this transition, it 
is helpful to consider the well-known mapping between 
phase transitions in c?-dimensional quantum systems and 
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FIG. 2: (a) The scattering processes transferring pairs of par- 
ticles between the two subbands. (b) The flow diagram for 
the renormalization group (JSJl. 



(d + l)-dimensional classical models [14[. In particu- 
lar, the phase transition in the one-dimensional quantum 
model ([3]) is equivalent to that in the two-dimensional 
classical Ising model [Hj]. In this mapping the gap A 
becomes the inverse correlation length r~ x of the Ising 
model, and the scaling r c oc |T — T c | _1 , well-known from 
the exact solutions, is equivalent to Eq. ([5]). The rela- 
tion between these phase transitions can be made more 
explicit by noticing that our Hamiltonian ((4]) essentially 
coincides with the transfer matrix [n| of the Ising model 
near the transition point. 

In the discussion leading to the result ([5]), the inter- 
actions in the quantum wire were assumed to be very 
strong. To explore the fate of the gap as the interac- 
tion strength is reduced, we now turn to the case of weak 
interactions. In this case the transition to the quasi- 
one-dimensional state occurs when electrons start popu- 
lating the second subband of transverse quantization in 
the wire. In this regime one can neglect the presence of 
the other subbands and present electron wavefunctions 
as ip(x,y) = ipi(x)xi(y) + ^2(x)x2(y), where xiAv) are 
the first and second eigenstates in the confining potential. 
Weak interactions between electrons lead to coupling of 
the two subbands. The low energy properties of the sys- 
tem are described by four interaction constants. The first 
three constants, g\, <?2, and g x , correspond to density- 
density interactions in the first subband, the second one, 
and between the two subbands, respectively. The fourth 
coupling constant g t accounts for the possibility of trans- 
fer of pairs of electrons from one subband to the other, 
Fig. Ha). 

It is well known that in one-dimensional systems 
density-density couplings renormalize the velocities of the 
acoustic low-energy excitations, but do not lead to the 
emergence of spectral gaps On the other hand, the 
coupling g t creates and destroys pairs of electrons in each 
subband and thus can, in principle, lead to a BCS-likc 
gap in the spectrum. Since the coupling constants in 
weakly interacting one-dimensional electron systems ac- 
quire logarithmic renormalizations at low energies 17j . 
the existence of a gap is determined by the scaling of gt- 

The renormalization group equations for the four cou- 



pling constants have been derived in Ref. [l8|. As the 
bandwidth of the problem is scaled from Do down to 
D, the renormalization of the coupling constants can be 
found by solving the system of two coupled equations 

y't = yyt, y' = yt- (6) 

Here the derivatives are with respect to £ = \n(Do/D), 
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vfi and vf2 are the Fermi velocities in the two subbands. 

The renormalization group flow corresponding to the 
equations ([6|) is shown in Fig. [2](b) . To find the ini- 
tial values of yt and y we compute the coupling con- 
stants in first order in the interaction strength. Assum- 
ing that the Coulomb interactions between electrons are 
screened by a gate at distance d, in the limit of low 
electron density in the second subband we find with 
logarithmic accuracy g^ — g x ^ — 2e 2 hiikpid) and 
g^' ~ e 2 (kp2d) 2 ln(l /kp 2 d). It is important to note that 
02 / V F2 vanishes when approaching the transition. This 
is a consequence of the Pauli principle. When the average 
distance between electrons is large, interactions between 
them are effectively local. Then, as identical fermions 
never occupy the same place, electrons essentially do not 
interact. From these estimates we conclude that j/ ' is 
positive, j/ -* ~ (3e 2 /7rfi.fFi) ln(fcFi^), and according to 
the flow diagram Fig. [DJb), the interaction constant gt 
scales to infinity. Consequently, the system develops a 
spectral gap. To find its value, we estimate g[ ^ near the 



transition as g t ~ e fc^/^Fi and obtain 
A oc (/j. - ii c ) a . 



(7) 



Thus the gap in the spectrum of transverse excitations of 
the wire exists not only when the interactions are strong, 
but also when they are weak. Unlike the case of strong 
interactions ([5]), at weak coupling the power-law depen- 
dence © has a very large exponent a = (4j/^ ^) _1 3> 1. 

To gain insight into the evolution of the transition 
between the two limiting cases, we derive the effective 
Hamiltonian of the system at intermediate interactions. 
Since the interactions between electrons in the lower sub- 
band are no longer weak and only their properties near 
the Fermi level are important, it is convenient to use the 
bosonization approach (l7j . On the other hand, as we dis- 
cussed, near the transition interactions between electrons 
in the second subband are negligible, 172 — * 0. Further- 
more, the curvature of their spectrum is important in this 
regime, so the description in terms of fcrmionic operators 
is more appropriate. The non- vanishing density-density 
interactions can still be described by the constants g\ and 
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g x , although their values may no longer be computed in 
first-order perturbation theory. Under these conditions 
the Hamiltonian has the form 

+7t J dx {[(dV0V> - tlidi>]e 2lK6{x) + h.c.j . (8) 

Here the bosonic fields <j)(x) and 9{x) describe the density 
excitations in the first subband, K = (1 + gi/nhvpi)^ 1 ^ 2 
is the respective Luttinger liquid parameter, tjj is the elec- 
tron destruction operator for the second subband, the 
constant jt ~ e 2 j and K = 1 — K 2 g x /i:hvpi. In deriv- 
ing Eq. © we performed a unitary transformation [2(| 
which removed the density-density coupling between the 
two subbands and changed the phase factor e 2ie in the 
last term to e 2lK - e . 

The Hamiltonian (JSj) interpolates between the limits 
of weak and strong interactions. In the weak coupling 
limit k — > 1, and a simple scaling analysis recovers the 
renormalization group equations ([6]). At strong interac- 
tions the problem of evaluating the coupling constants g\ 
and g x is non-trivial, but in the regime e 2 /hvpi <C 1 <C 
(e 2 /Hvpi) hi(kpid) one can still use our earlier estimates 
gi = g x = 2e 2 ln(fcpi<i) and conclude that k — > 0. Inter- 
estingly, in this case the bosonic and fermionic parts of 
the Hamiltonian ((SJ) decouple, with the latter becoming 
equivalent to Eq. (J4j) . 

One can use the Hamiltonian JHJ) to discuss the evolu- 
tion of the gap with varying interaction strength. In the 
limit of strong interactions, when k = 0, the magnitude 
of the pairing term scales to zero near the transition as 
kp^, i.e., slower than the Fermi energy Ep 2 oc k F2 - In 
this regime, the gap equals the Fermi energy, A = Ep 2 , 
cf. Eq. At strong but finite interactions the pairing 
term suffers additional power-law suppression at kp 2 — > 
because of the factor e 2tKe . However, as long as it scales 
slower than the Fermi energy, the magnitude of the gap 
remains A = Ep 2 . At weaker interactions, when k ex- 
ceeds a certain critical value, the pairing term scales to 
zero faster than the Fermi energy. In this regime the 
gap develops in a small vicinity of the Fermi points, and 
its dependence on chemical potential is given by a non- 
universal power law ([7]) with exponent a > 1. A more 
detailed theory of the transition at intermediate interac- 
tion strengths will be reported elsewhere [2lj |. 

Our results are summarized in the phase diagram 
Fig. [T] The electron system in a quantum wire re- 
mains one-dimensional and has a single acoustic excita- 
tion branch until the chemical potential reaches a certain 
critical value \x c . At the critical point there is a second 
gapless mode, and the system can no longer be viewed 
as one-dimensional. At /1 > [i c the second mode devel- 
ops a gap A oc (/i — /i c ) Q with exponent a = 1 at strong 



interactions, but very large a at weak coupling. At weak 
interactions, as the chemical potential is increased fur- 
ther, the residual interactions g 2 grow, and the gap dis- 
appears. This happens when the electron density in the 
second subband is still very small, of order (kpd 2 ) -1 , 
long before the third transverse mode appears. We see 
no physical reason for the gap to disappear at strong cou- 
pling. In experiments, the presence of a gap will affect 
the temperature dependence of the conductance which is 
expected to show activated behavior even above the tran- 
sition into the quasi-one-dimensional state. The doubling 
of the zero-temperature conductance (from e 2 /h to 2e 2 /h 
for spinless electrons) occurs at the upper line in Fig. [T] 
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